SHARP ESTIMATES OF THE JACOBI HEAT KERNEL 



ADAM NOWAK AND PETER SJOGREN 

Abstract. The heat kernel associated with the setting of the classical Jacobi polynomials is 
defined by an oscillatory sum which cannot be computed explicitly, in contrast to the situation 
for the two other classical systems of orthogonal polynomials. We deduce sharp estimates giving 
the order of magnitude of this kernel, for type parameters a, f3 > —1/2. As an application of 
the upper bound obtained, we show that the maximal operator of the multi-dimensional Jacobi 
heat semigroup satisfies a weak type (1,1) inequality. We also obtain sharp estimates of the 
Poisson- Jacobi kernel. 



1. Introduction 

Let Pn^i n = 0, 1, 2, . . ., be the classical Jacobi polynomials with type parameters a, f3 > — 1, 
as defined in Szego's monograph |23j. The Jacobi heat kernel is given by 

(1) G^(s,y)=^exp(-tn(n + a + £ + l)) " [X) a J [V) , x,y€[-l,l], t > 0, 

n=0 ' "n' 

where hn^ = /_i [-Pn'^ \ x )] 2 (1 — x) a (l + x)@dx are proper normalizing constants. The numbers 
n(n + a + f3 + 1) here are the eigenvalues of the Jacobi differential operator, and it is well known 
that this kernel gives the solution of the initial-value problem for the Jacobi heat equation, as 
explained below. 

Our main result reads as follows. 

Theorem A. Assume that a,/3> —1/2. Given any T > 0, there exist positive constants C, c\ 
and C2, depending only on a, ft and T, such that 

1 T + , a l-a-V2r. . / flW m -0-1/2 1 ( (0 - <f) 2 \ 
-{t + 6ip\ [t + (ir - 6){ir - tp)\ -^=exp^-ci I 

< G"'^(cos 9, cos (p) 

<C[t + 9<p] — 1/2 [t+{*- 9)(n - v) ] -"- 1/2 1= exp (-cJ^^D) , 

for 9, <p G [0, 7r] and < t < T. Moreover, 

C- 1 <G^(x,y)<C, 

for x, y G [—1, 1] and t > T, and Gf (x, y) — > I/Ziq'' 3 as t — > oo, uniformly in x,y € [—1, 1]. 

Thus we obtain a qualitatively sharp description of the behavior of Gf (x, y). The restriction 
on a and (3 is imposed by the methods used. Nevertheless, it is natural to conjecture that the 
same bounds hold for all a,/3 > —1. 
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The multi-dimensional Jacobi heat kernel is a tensor product of one-dimensional kernels, and 
Theorem |A] provides similar bounds also in the multi-dimensional setting. As an application, 
we prove that the maximal operator of the multi-dimensional Jacobi semigroup satisfies a weak 
type (1,1) estimate, see Theorem 15. 11 This complements analogous results in the Hermite and 
Laguerre polynomial settings, which were obtained in dimension one by Muckenhoupt |14] , and 
in arbitrary finite dimension by the second author [21] and Dinger [8], respectively. For the 
Laguerre case, see also the authors' paper [XT] . 

The heat kernels associated with the two other families of classical orthogonal polynomials 
have been known explicitly for a long time. Already in 1866, Mehler |13j established the formula 

H n (x)H n {y) ^ n _ 1 ^ / 2xyr - (x 2 + y 2 )r 2 

71=0 

which makes it possible to sum the heat kernel related to the Hermite polynomials H n . In the 
case of the Laguerre polynomials L", the relevant bilinear generating function is the Hille-Hardy 
formula 

t ,. „ t:mM = ^- exp ( - M+yt) M -^, a ( M , 

1 (n + a + 1) 1 — r \ 1 — r J \ 1 — r J 

where \r\ < 1, a > — 1 and I a is the modified Bessel function of the first kind. This identity 
was found in 1926 by Hille [TT] and independently rediscovered later by Hardy [ 10] , see [25]. An 
analogue of these formulas in the Jacobi setting is Bailey's formula 

^ (cos fl)^ (cosy) r(q + /3 + 2) 1-r 

h a,p r 2 °+/ 3 + 1 r(Q + i)r(/3 + 1) (1 + r)-+/3+ 2 

/o+^+2 a + /3 + 3 / 2 sin | sin f \ 2 / 2 cos f cos f \ 2 \ 

XF \ 2 ' 2 > a + 1 >P + 1 > [ r l/2 + r -l/2 ) > ^1/2 + r -l/2 ) J' 

where |r| < 1, a,/3 > —1 and is Appel's hyper geometric function of two variables. This 
generating function was first stated in 1935 without proof in Bailey's tract [lj. The proof is 
a straightforward consequence of Watson's formula for hyper geometric functions [24] and was 
published slightly later [2]. 

However, in contrast with the Hermite and Laguerre cases, Bailey's formula does not enable 
one to compute the Jacobi heat kernel. This is because the eigenvalues n(n + a + /3 + l) occurring 
in the defining series are not linear in n. It is known that in the four simple special cases a, /3 = 
±1/2, the kernel G?' p (x,y) can be written by means of non-oscillating series. The argument 
is based on the periodized Gauss- Weierstrass kernel and simple initial-value problems for the 
classical heat equation in an interval. No further elementary representation for G^ =1//2 ' ±1//2 (x, y) 
seems to be possible. This indicates that there is little hope of deriving a closed formula for 
the Jacobi heat kernel for general a and (3 similar to those of Mehler, Hille and Hardy, and 
Bailey. A natural and desirable substitute for an exact expression are therefore the estimates in 
Theorem [A] 

The Jacobi polynomials cover as special cases several other classical families of orthogonal 
polynomials, including Chebyshev, Legendre and Gegenbauer (also called ultraspherical) poly- 
nomials. Special instances of the Jacobi heat kernel exist at least implicitly in the literature 
since the 19th century, and the question of describing its behavior was an open problem, even 
though perhaps never stated explicitly in written form. Additional motivation comes from the 
fact that Gf (x,y) is also the transition probability density for the Jacobi diffusion process, 



\r\ < 1, 
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which has important applications in stochastic modeling in physics, economics and genetics; see 
[12] and references given there. According to our knowledge, the behavior of Gf (x,y) has not 
been investigated before, except for its positivity. Bochner [3] proved that the ultraspherical 
heat kernel is non-negative. Strict positivity in the general Jacobi case was shown by Karlin and 
McGregor [12]. Some later results on the positivity can be found in Gasper [9] and Bochner [3]. 

We also take the opportunity, see the Appendix, to describe the behavior of the Poisson-Jacobi 
kernel, which is essentially the sum occurring in Bailey's formula. However, the representation 
in terms of Appel's function F± does not seem to be very useful for this purpose. Instead we 
employ a double integral representation that was derived recently by the authors [20] from a 
product formula due to Dijksma and Koornwinder [7]. 

For short times t, a direct analytic treatment of the heavily oscillating series defining Gf (x, y) 
is practically impossible. Therefore, we develop a method combining several ingredients. These 
are, among others, the already mentioned product formula of Dijksma and Koornwinder and a 
resulting reduction formula, transference of heat kernel estimates from a sphere, a comparison 
principle relating heat kernels for different type parameters, the semigroup property and, finally, 
a rough estimate of the series defining Gf'^(x,y). 

The paper is organized as follows. In Section [2] we introduce three related Jacobi settings 
appearing in the literature, and explain how the associated heat kernels are connected. Section 
[3] contains the auxiliary results that form the main tools of our proof of Theorem [A] Several 
of them are interesting in their own right. Section [4] is devoted to the proof of Theorem [A] 
Maximal operators of multi-dimensional Jacobi semigroups are treated in Section [5j by means 
of Theorem [Al we prove weak type (1, 1) estimates for these operators. Finally, in the Appendix 
we prove sharp estimates for the Poisson-Jacobi kernel. 

Throughout the paper we use standard notation. The letter C will stand for many different 
positive constants independent of significant quantities. When writing estimates, we use the 
notation X < Y to indicate that X < CY. We write X ~ Y when simultaneously X < Y and 
Y < X. Tracing the proof of Theorem it is easy to verify that the constants in the statement 
come out as claimed. 



2. Preliminaries 

Given a,/3> —1, the one-dimensional Jacobi polynomials of type a,/3 are defined for n £ N 
and — 1 < x < 1 by the Rodrigues formula (cf. [23], (4.3.1)]) 



n v i 2 n n\ y 1 K ' dx n 



\3+n 



Note that each P" is a polynomial of degree n. 

We will consider three closely related settings of orthogonal systems based on Jacobi poly- 
nomials. All of them have deep roots in the existing literature. Below we briefly introduce 
each setting, for the sake of simplicity in dimension one. Multi-dimensional analogues arise in a 
standard way as tensor products of one-dimensional systems. 

Pure polynomial setting. In this case the relevant system {P" : n > 0} is formed directly 
by Jacobi polynomials. This system is an orthogonal basis in L 2 (dQ a ^), where g a ^ is the 
beta-type measure given by 

dQaA x ) = (1 - xTi 1 + x f dx 
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in the interval [—1,1]. Each Pn is an eigenfunction of the Jacobi differential operator 

jocfi = _ (1 _ x ^&_ _ ^ _ a _ (a + p + 2)x i d 
ax* dx 

more precisely 

J a >Pp^ = n(n + a + /3 + l)P%' 13 , n > 0. 

The operator J a ^ is symmetric and nonnegative in L 2 {dQ a ^) on the domain C 2 (— 1,1), and 
has a natural self-adjoint extension whose spectral resolution is given by the Pn , see [18] for 
details. The semigroup T" ,/3 = exp(— iJ Q '^) is a symmetric diffusion semigroup in the sense of 
[22} Chapter 3]; in particular, Tj*'^l = 1. It is also the transition semigroup for the Jacobi 
diffusion process, which already received attention; cf. [12] and references there. Some aspects 
of harmonic analysis in the multi-dimensional Jacobi pure polynomial setting were investigated 
by the authors in [18] . 

The integral representation of T™ , valid for / G L l {dg a ^), is 

T t a '?f(x) = j 1 ^(x,y)f(y)dg a ^y), 

where the Jacobi heat kernel is given by the oscillating series ([I]). The normalizing constants 
h a n 3 := \\Pn^\\ 2 L 2 (dea ip) in © are given by (cf. [23 (4.3.3)]) 

(2) h *j, = 2 a +/ 3 + 1 r(n + q + l)r(n + /3 + l) 

{ ' n (2n + a + (3 + 1)I> + a + p + l)T(n + 1) ' 

where for n = and a + j3 = — 1 the product (2n + a + /3 + l)r(n + a + /3 + 1) must be replaced 
by r(a + /3 + 2). Notice that l//in'^ — n, n > 1. As already mentioned, G^'^{x,y) is strictly 
positive for x,y G [—1,1] and £ > 0. Further, G" ,/3 (x,y) is continuous (and even smooth) in 
(t,x,y) G (0, oo) x [— 1,1] 2 , as well as in (a,/3) G (— l,oo) 2 . This can be verified by analyzing 
the defining series, using the corresponding continuity properties of Jacobi polynomials and the 
bound 

\P^(x)\<(n + lf, n>0, xG[-l,l], 

justified in Section [3l By [TBI Proposition 3.3], f{x) is for any / G L l {dg a ^) a C°° function 
of (t,x) G (0, oo) x (—1, 1) satisfying the heat equation 

(d t + J a >P)T t a >Pf(x) = 0, xG(-l,l), t>0. 

Finally, we note that (see [IFJ p. 347]) for / G C[— 1, 1] 

lim T^/(x)=/(x), xG [-1,1], 

i->0+ 

and the convergence is uniform in x. All these facts will be used in the sequel. 



Trigonometric polynomial setting. This framework emerges if one applies the natural and 
convenient trigonometric parametrization x = cos 9, 9 G [0,7r], to the Jacobi polynomials. We 
consider the normalized trigonometric polynomials 

V^{9) = 2( Q +^+ 1 )/ 2 (/CT 1/2p n* /3 (cosfl). 

The system {Vn ,/3 : n > 0} is orthonormal and complete in L 2 (dfi a ^), where 

/ 6»\2q+1/ 6\ 2/3+1 
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in [0,7r]. Each Vn^ is an eigenfunction of the differential operator 

r a,i3 _ d<2 a - (3 + (a + (3 + l)cos8 d fa + f3 + l\ 2 ^ 



d9 2 sin0 dO \ 2 

indeed 



J a ' p ^' p = (n + j K' P , n > 0. 

The operator jT a '^ has a natural self-adjoint extension whose spectral resolution is given by 
the Vn^ , see [20] for details. The semigroup T t a ^ = exp(—tj' a '^) has the integral representation 

7^/(0) = ¥>)/(¥>) d/W^). 

JO 

valid for / 6 L 1 (d/i Q , i( g), with the heat kernel defined by 

oo 
n=0 

Note that J7 a '^ is obtained by transforming J"'' 3 according to the change of variable x = cos 
and introducing the zero order term. The latter modification leads to eigenvalues which are 
squares, and therefore the oscillating series defining the one-dimensional Poisson-Jacobi kernel 
in this setting, 

oo 
n=0 

can be represented in a more convenient way; in particular, Bailey's formula applies. On the 
other hand, the semigroup T^'^ is submarkovian, but not Markovian in general. 

Several fundamental harmonic analysis operators related to the Jacobi trigonometric poly- 
nomial setting were studied recently by the authors [20J. The ultraspherical case was widely 
investigated from a slightly different perspective in the seminal paper of Muckenhoupt and Stein 
|16j . which in 1965 initiated the development in harmonic analysis known as harmonic analy- 
sis of orthogonal expansions. In both cases, the analysis was based on the one-dimensional 
Poisson-Jacobi kernel; see [201 for further facts and references. 



Trigonometric 'function' setting. This context originates naturally in connection with trans- 
plantation problems for Jacobi expansions (see for instance [151 15] and references there) and is 
derived from the previous setting by modifying the Jacobi trigonometric polynomials so as to 
make the resulting system orthogonal with respect to Lebesgue measure d6 in [0,7r]. More 
precisely, we consider the functions 

C /3 W = (sin-) (cos-) V^(9), n>0. 

Then the system {4>n'^ '■ n > 0} is an orthonormal basis in L 2 (d9). The associated differential 
operator is, cf. [231 Section 4.24], 

W> - _iL + («-l/2)(« + l/2) 09 - 1/2)09 + 1/2) 



and we have 



d9 2 4 sin 2 f 4 cos 2 f 



rV=fn+^±|±^)V, n>0. 
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The semigroup T"'^ = exp(— iJP'' 3 ), generated by the natural self-adjoint extension of J a 'P , 
has the integral representation, valid for / G L 2 (d6), 



T^f{9)= / G^(e,<p)f(<p)d<p, 
Jo 

where 

oo 



n=0 



Note that T" ,/3 is not defined on all U\dO), 1< p < oo, if a < -1/2 or < -1/2. 
The Poisson-Jacobi kernel in the Jacobi trigonometric 'function' setting is defined by 



n=0 



Observe that there is a simple analytic connection between the heat kernels in the three Jacobi 
frameworks. In fact, we have 



(3) = 2 Q+/3+i e H 2 M sin -sin -J (cos -cos -J G"' p (cos (9, cos tp). 

V 2 2/ V 2 2/ 

Similarly, for the Jacobi-Poisson kernels, 

(^.W = (^sm-sin-J (^cos-cos-J H t (0,<p). 
Thus kernel estimates can be translated between the Jacobi settings. 

3. Preparatory results 

In this section, we prove several results that will be important ingredients of the proof of 
Theorem 1X1 They are of independent interest, and so some of them are stated in slightly larger 
generality than actually needed for our present purposes. 

It is convenient to introduce a compact notation for objects related to the ultraspherical 
setting, i.e., when the Jacobi parameters are equal, say a = /3 = A. In such cases, the sub- or 
superscript A, A will be shortened to A; for instance 

P n '■= P n' X i QX '■= QX,Xi '■= h n X - 

Notice that this convention differs somewhat from the standard notation for the classical ultra- 
spherical (Gegenbauer) polynomials C„, cf. [231 Section 4.7]. In fact we have, see [231 (4-7.1)], 

U) c A fx) _ r(A + l/2) r(n + 2A) pX - 1/2() . 

(4) Cn{x) - r(2A) r(n + A + l/2) P " (X) ' A> -V2, A ^ 0. 

Reduction formula. The following product formula for Jacobi polynomials was derived by 
Dijksma and Koornwinder [7]: 

W»M oJwPh r(a + /3 + l)r(n + a + l)r(n + /3 + l) 



im\T(n + a + (3 + l)T(a + l/2)r(/3 + 1/2) 
11 C^ P+1 {ust + vs/l^s/l^){\ - u 2 Y- l ' 2 (l - v 2 )^ l ' 2 dudv. 



-l J-i 
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This formula is valid for a, f3 > —1/2. We shall write it in a more suitable form which, by a 
limiting argument, will be valid for all a, /3 > —1/2. 

Let Tl a be the probability measure on the interval [—1,1] defined for a > —1/2 by 

du a ( U )= if + 1} (i - u r- i/2 du. 

V7iT(a + 1/2) 



In the limit case a = —1/2, we put 



n_i/ a = ^(A-i + tfi), 



where 5±i denotes a point mass at ±1. Note that n_ 1 / 2 is the weak limit of H a as a — > —1/2 

| and t = sin % 



We now rewrite the above product formula with s = sin § and i = sin % , using the relation 



between and K ,/3 , the fact that (cf. [231 (4.1.1)]) 



A . 1 ._ T(n + A + 1) 



r(n+i)r(A + i) 

and the expression ([2]) for hn'^ ■ After some computations, one finds that 

_y/*T{a + + 3/2) tr . vm . P 2 ^ +1/ ^nsin|sin|+ W cos|cos|)P 2 a n +/?+1/2 (l) 



J I dn a {u)dnp(v) 



r(a + i)ro? + i) J J aK ' pw C m/2 

This formula holds for all a,0 > —1/2, since Jacobi polynomials are continuous functions of 
their type parameters (see [231 (4-21.2)]). 

Multiplying both sides above by exp(— t(n + Q+ 2 +1 ) 2 ) and summing over n > we get 



oo p a+(3+l/2 ( . ■ <£_> £\ p a+/3+l/2 n , 

-tn(n+Q+/3+l) f2n ^ Sm 2 Sln 2 + V COS 2 COS 2 > r 2n ' 1 ' 

7 ' a+/J+l/2 



n=0 "2n 



Writing tn(n + a + /3 + 1) = |2n[2n + (a + /3 + 1/2) + (a + /3 + 1/2) + 1] and taking into 
account that ultraspherical polynomials of even (odd) orders are even (odd) functions (cf. [23[ 
(4.7.4)]), we see that the last series represents the even part with respect to the first variable 
of the ultraspherical heat kernel with the second variable fixed at the endpoint 1. Since for 
symmetry reasons the corresponding odd part gives no contribution to the integral, we end up 
with the following reduction formula. 

Theorem 3.1. Let a,/3> -1/2. Then, for all 6,(p G [0,tt] and t > 0, 

Gt ,/3 (cos 6, cos <p) = C at p f f G"f/ +1/2 (u sin ^ sin ~ + vcos ^cos^, l) dIl a (u)dILp(v), 



*/4 V 2 2 2 2 

with the constant C a $ = ^T(a + (3 + 3/2)/(2 Q +^ +1 r(a + l)T(f3 + 1)) . □ 

Thus we have expressed the general Jacobi heat kernel in terms of the ultraspherical one. The 
relation between the type parameters here will be essential for our arguments to prove the heat 
kernel estimates. 
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Connection with the spherical heat kernel. We now consider the Jacobi polynomial setting 
in the half-integer ultraspherical case, that is when 

N 

a = /3 = --l, N = l,2,.... 

It is well known that this situation is closely connected with expansions in spherical harmonics 
on the Euclidean unit sphere of dimension N, see for instance |12|, Section III]. In particular, 
there exists a relation between the heat kernels in the two settings, which we indicate below. 

For N > 1, let S N be the unit sphere in M Ar+1 and denote by ajq the standard (non-normalized) 
area measure on S . The Laplace-Beltrami operator A^v on S N is symmetric and nonnegative in 
C°°(S N ) C L 2 (da^). The classical system of spherical harmonics on S N is an orthogonal basis 
in L 2 (cI(Tn) of eigenfunctions of A^v- The spherical heat semigroup J7 t = exp(tATv), generated 
by the self-adjoint extension of Ajv, has an integral representation 

U t N f(0= I K?{j-,ri)f{n)d*T N {rj), $eS N , t > 0, 
Js N 

for L 2 (dcJiv). The spherical heat kernel K% can be expressed explicitly as an oscillatory 

series of spherical harmonics, see \12\ p. 176]. By general theory (cf. [6, Theorem 5.2.1]), 
Kt (£j v) is a strictly positive and continuous (even smooth) function of (t, £, rf) 6 (0, oo) x S x 
S N . 

It is well known that the zonal case in the context of A^y and expansions with respect to 
spherical harmonics reduces to the ultraspherical setting in the interval [—1,1] with the type 
parameter A = N/2 — 1. Indeed, let F be a zonal function on S , say F = f o ip f where 

is the zonal projection onto the diameter of S N determined by the first coordinate axis. Then the 
expansion of F in spherical harmonics reduces to the expansion of / in ultraspherical polynomials 
of type A = N/ 2 — 1. The associated heat semigroups are related in a similar way, as stated 
below. 

Lemma 3.2. Assume that N > 1 and A = N/2 — 1. Then we have for f G L 2 (dg\) 

{T t x f)om = U t N (fo^)(0, £eS N , t>0. 

Proof. Observe that the semigroups considered consist of L 2 -bounded linear operators, which 
for each t > map I? into the subspaces of continuous functions on [— 1 , 1] and S N , respectively. 
Moreover, / G L 2 (dg\) if and only if / o ip e L 2 (do~N). Therefore, since linear combinations of 
ultraspherical polynomials are dense in L 2 (dg\), we may assume that / = P fc A for some k. 

The identity we must prove has roots in the fact that the ultraspherical operator J A is essen- 
tially the zonal part of the Laplace-Beltrami operator A^r on S N . Indeed, writing the differential 
operator Ajy in hyperspherical coordinates on 5^ (see jT2l p. 175]) one easily verifies that 

(j A P fc A )o^ = -A w (P fc A o^). 

Since P£ is an eigenfunction of J A , we see that P^ o tp is an eigenfunction of A^r, with the same 
eigenvalue. For smooth functions, the self-adjoint extension of A^v coincides with the differential 
operator, and we conclude that 

(T A P A ) oiKfl = e-WA+DpA Q m = tyiV (p A ^ (e)) f £ S N > 

as desired. □ 
We now establish a connection between the ultraspherical and spherical heat kernels. 
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Theorem 3.3. Assume that N > 1 and A = N/2 - 1. T/ien 

G^y) = J i (e, (y, CVl " y 2 )) <^jv-i(C), *,yG[-l,l], t > 0, 
where £ G is arbitrary. 

Proof. Let / be a polynomial on [—1,1]. By Lemma 13.21 



W(£),I/)/fo)(l-I/T<fo= / iff (C,r?)/ oip( v )do- N ( V ). 

Js N 

To treat the last integral, we introduce zonal coordinates on S N , 

¥: [-1,1] x S^ 1 — ► S N , *(y,C) = (^Cv 7 ! 3 ^)- 
Then for reasonable functions -F 

/ F(0da N (0= ^ f FoM ) (y,()da N ^(0(l-y 2 ) N/2 - 1 dy. 



Therefore, 



1 <#ty(0,i/)/(i/)(i - y 2 ) x dy = f 1 f ^ K?(t, {yxV^y 2 )) ^-i(C) /(v)(i - y 2 ) A */- 

Since polynomials are dense in C[— 1, 1] and the kernels in question are continuous functions of 
their arguments, it follows that 

JSN-l v / 

for t > 0, y G [-1, 1] and £ G S N . □ 



The case y = 1 and £ = (x, Vl-x 2 , 0, . . . , 0) of Theor em 13.31 reveals a particularly simple 
relation between the ultraspherical and spherical heat kernels. 

Corollary 3.4. Let A = N/2 - 1 for some N G {1, 2, . . .}. Then 

Gf(x, 1) = o-N^iS"- 1 ) K t N ((x, Vl - x 2 , 0, . . . , 0) , (1, 0, . . . , 0)) , xe [-1,1], t > 0. 

This expression for the ultraspherical heat kernel in terms of the spherical one will allow us 
to transfer qualitatively sharp heat kernel bounds on spheres to the ultraspherical setting with 
half- integer type index. On the other hand, it is interesting to observe that the spherical heat 
kernel on S N is completely determined by G^(x,l) for A = N/2 — 1. This is a consequence 
of Corollary 13.41 and the fact that K^(^,rj) depends on £ and r\ only through their spherical 
distance. 

Comparison principle. Given e, S > 0, define 

^s(x) = (1 - xY' 2 {l + xf 2 , x G [-1, 1], 

with the convention that (l±x) = 1 for x = =pl. This is the square root of the Radon-Nikodym 
derivative dg a+€> p+s I dg a ,p ■ Using a parabolic PDE technique, we shall prove the following result 
comparing Jacobi heat kernels with different type parameters. 

Theorem 3.5. Let a,/3 > — 1. Given e,5 > and a > — e/2, f3 > —5/2, we have 



(5) * £ , a (s)<My) Gr ,P+ >.2/) <exp(— (a + /3 + l + — )t) G? p (x,y) 

for all x, y G [—1,1] and t > 0. 
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Translating this estimate to the other two Jacobi settings, we get 
Corollary 3.6. Let a, /3,e,8 be as in Theorem \3.5[ Then 

( sin \ sin |) ( cos \ cos |) G? +t ' p+s {0, < £f '^(A, tp) 

and 

for all 9, if £ [0, tt] and t > 0, with the natural convention for boundary values in the second 
inequality. 

The relation in the Jacobi trigonometric 'function' setting is particularly nice since it shows 
that the heat kernel is decreasing as a function of each of the type parameters a > and /3 > 0. 
Note that, by subordination, the estimates of Corollary 13.61 carry over to the corresponding 
Poisson kernels. 

Proof of Theorem 13.51 Since the Jacobi heat kernel is a continuous function of its type param- 
eters a,/3 > —1, we may assume that e, 5 > 0. We first rewrite ([5]) in integrated form. By 
integrating against f(y) dg a ^(y), we see that ([5]) implies 

(6) ^ s (x)T t a+ ^ + \f/^,s)(x) < e^W+l+^t T ^f{x) 

for suitable functions / > 0. Conversely, if ([6]) holds for all nonnegative / £ C£°(— 1,1), then 
© will follow, since Gf^(x,y) is continuous in (x,y) 6 [—1, l] 2 . We shall thus prove the lemma 
by verifying © for x £ (—1,1) and t > and any < / G C£°(— 1,1) not identically 0. Our 
reasoning will rely on a generalization of the minimum principle method used to prove |18} 
Lemma 3.4]. 

Denote by u = u(t, x) the left-hand side of © and let 

V = V (t, X) = e tV e ^(a + fi+l+^)t T «,P {f + v){x) 

for some fixed r/ > 0. Since / is smooth, the functions u and v have continuous extensions to 
[0, 00) x (—1, 1). Our task will be done once we show that 

u(t,x) < v(t,x), x £ (—1,1), 

for all t > and any i] > 0. Let 

T = sup{t' > : u(t,x) < v{t,x) for (t,x) G [0,i') x (-1,1)}. 

Clearly, u(0,x) = f{x) < f(x) + rj = v(0,x) for x £ (—1,1). Moreover, u(t,x) < v(t,x) for all 
t > provided that \x\ is sufficiently close to 1; this is because u(t, x) < C<& e ^{x) and v(t, x) > r] 
for t > 0, x £ (—1, 1). Hence for t small enough u(t,x) < v(t,x), x £ (—1, 1), which means that 
T > 0. 

Suppose that T is finite. We shall then derive a contradiction, which will end the reasoning. 
Observe that u(T, x) < v(T, x) for all x £ (—1, 1) and u(T, xo) = v(T, xq) for some xo £ (—1, 1). 
We claim that 

(7) d t (v(t, x) - u(t, z)) I {t)X)={T)Xo) > 0. 

This would imply that v(t, xq,) — u(t, xq) < for t slightly less than T, a contradiction. 

To prove this claim, we compute the derivative in ([7]). With the aid of the heat equation we 
get 

d t (v(t,x) - u(t,x)) = 
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6 + 6 (a + + 1 + ^) + 7/J v(t, x) - J a ^v(t, x) + <S> e , 8 (x)J a+ ^ +s (u(t, x)/* e 4x)) . 

Then using the definition of J a ^ and the fact that v — u = d x (v — u) = at the point (T,xq), 
we find after somewhat lengthy computations that the left-hand side in (|7|) is equal to 



(l- x l)d 2 x (v-u)(T,x Q ) + 



e(a + e/2) | 5(0 + 5/2) 

1 — X 1 + ^0 



u(T, Xq) + 7]U(T, Xq). 



The first term above is nonnegative, since the function x \— > v(T, x)—u(T, x) has a local minimum 
at x = xq. The factor in the square bracket is obviously nonnegative by the assumptions 
on e, 5, a, /3. Finally, u(T, xq) is strictly positive by the corresponding property of the kernel 
involved. The claim follows. □ 

We remark that when either a < —e/2 and e > or f3 < —5/2 and 5 > 0, the estimate of 
Theorem 13.51 (and thus also the estimates of Corollary I3.6[) does not hold. This can be shown 
by means of a counterexample very similar to that of [181 Remark 3.6]. 

Rough estimate. We now employ absolute value estimates of Jacobi polynomials to obtain a 
rough short time bound for the Jacobi heat kernel in terms of t only. This method, of course, 
distinguishes no subtle effects coming from oscillations. Therefore, the resulting estimate is far 
from sharp. More accurate upper bounds for the Jacobi heat kernel, involving also dependence 
on x and y, can be found by means of a more detailed analysis and the estimates for Jacobi 
polynomials contained in [23|. Theorem 7.32.2]; see also [231 (7. 32. 6), (7. 32. 7)]. 

Theorem 3.7. Let a,/3 > — 1 and T > be fixed. Then 

C/'\x,y)<t- c \ x^G [-1,1], 0<t<T, 

where the constant Co depends only on a and (3. 

In the proof we will use the following bound for Jacobi polynomials (see [23, (7.32.2)]) 

(8) \P^(x)\<n\ n>l, a; €[-1,1], 

where 7 = max{a,/3, —1/2}. 

Proof of Theorem\3.1\ Recall that 1/hn — n f° r n> 1. Thus 



G^(x,y) <l + J2 e ~ tnin+a+P+1)n \ P n^^)\\ p ^(y)\ <l + J2e- tn n\P^(x)\\P^(y)\, 



n=2 n=2 

and (jSJ) implies 



00 1 00 

G?'?(x, y) < 1 + £ e-V+* < 1 + ^ J2 e~ tn/2 < t^~\ 

n=2 n=2 

The theorem is proved. □ 
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4. Proof of Theorem [X] 

The large time behavior of G^ ,/3 (x, y) stated in Theorem lAl is a consequence of the short time 
estimate, in the following way. Given T > 0, the short time bound implies that 

G«'P(z,y)~l, z,ye [-1,1]. 

By the semigroup property, for t > T and x,y £ [—1,1] one has 

Gt'^&v) = j^G^ T (x,z)G^{z,y)dg a ^(z). 

Since 

J ^ G?J T (x, z) dg a ,p{z) = T t %l(x) = 1, x E [-1, 1], t >T, 
we get the estimates 

G^(x,y)~l, x, ye [-1,1], t > T. 

The existence of the uniform limit as t — > oo follows by combining the oscillating series ([I]) with 
the estimate (|8|) for Jacobi polynomials. 

Thus it remains to prove the short time estimates, and we first introduce some further no- 
tation. A real number r will be called dyadic if r = n/2 k for some integers n,k. We will use 
the notation X ~~ Y exp(-cZ) to indicate that Y exp(-ciZ) < X < Yexp(— c%Z\ with posi- 
tive constants c\ and C2 independent of significant quantities. Thus the short time estimates of 
Theorem [A] can be written 

(9) G t a ^(cos0,cos^) 

/ 9 6» (p\-P-i/2 1 ( (9-vf 

~~ I i + sin - sin — J 1 1 + cos - cos — — = exp 



2 2/ V 2 2/ V* V * 

The quantity T > will be fixed for the rest of the proof of Theorem lAl For the sake of clarity, 
we divide the proof into several steps, as follows. 

1. Estimate G^(x,l) for half- integer A > —1/2 by transferring, via Theorem 13.31 known 
bounds for the spherical heat kernel. 

2. Starting from the bounds of Step 1, iterate the reduction formula (Theorem 13.1)) to 
estimate G£(x, 1) for all dyadic values of A > —1/2. 

3. Apply the reduction formula to the estimate of Step 2 to prove (J9j) when a,/3 > —1/2 
and the sum a + /3 is a dyadic number. 

4. Combine the estimate of Step 3 with the comparison principle (Theorem 13. 5p to obtain 
a weakened version of in the ultraspherical case for a = f3 = A > 0. 

5. Use the semigroup property and the rough estimate of Theorem 13.71 to eliminate the 
weakening in Step 4 and prove ([9]) fora = /3 = A>0. 

6. Use the estimate of Step 5 and iterate with the reduction formula to estimate G^(x, 1) 
for — 1/2<A<0. In the final stroke apply again the reduction formula to prove ([9]) for 
all a,p> -1/2. 

Step 1. We first invoke the well-known Gaussian bounds for the spherical heat kernel, see [BJ 
Theorems 5.5.6 and 5.6.1]. Let d(£,ri) = arccos(£, rj) be the spherical distance between £ and 
rj G S N . Given any 5 > 0, we have 

r N(e: x / 1 ( d(tV? \ * _^oN 



^ V (e,r/)<^exp^-^^), t^eS", 0<t<T, 
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and 



These estimates, together with Corollary 13.41 and the observation that d{£,rj) ~ |£ — rj\ for 
£,77 G S , lead to the following result. 

Lemma 4.1. Assume that A = -/V/2 — 1 for some N > 1. T/ien 

G£(x,l)~~p^exp(-c^Y xE[-l,l], 0<t<T, 

or equivalently, 

1 / 9 2 \ 

(10) G£ (cos 0,1) ~~ ^+T exp ( ~ c t)' 6 I ' 71 "]' 0<t<T. 

Notice that this estimate is a special case of ©. 

Step 2. We claim that the bounds of Lemma 14.11 hold for all dyadic values of A > —1/2. To 
verify this, it is enough to prove the following lemma, since one can then iterate. 

Lemma 4.2. Assume that the estimate (fTO)) holds for some A > —1/2. Then it holds also with 
A replaced by X = A/2 - 1/4. 

To prove this lemma, we need an auxiliary result. 

Lemma 4.3. Let v > -1/2. Then 

J exp(zs) <flL,(s) ~ (1 + z)- u - 1/2 e z , z>0. 

Proof. One can assume that z > 1, since the opposite case is trivial. We split the integral and 
observe that the integral taken over (0,1) is larger than that over (—1,0). Thus we need only 
consider 



/ e zs dU u (s) ~ f e z8 (l-sf- x / 2 ds. 
Jo Jo 



Here we make the two transformations t = 1 — s and r = zt, and get 

e^— 1 ' 2 [ Z e-W-^dr. 
Jo 

The conclusion follows. □ 
Proof of Lemma \4-2\ Suppose that (fl~0j) holds for some A > —1/2. By Theorem 13. 1| 

G t y (cos0,l) = Cn A ,([-l,l]) J G^ /4 (vcos^,l)dUy(v) 

1 f f 1 — V cos - \ 
(11) ~~ J exp I - c 2 -j dUy(v). 

Applying now Lemma [4.31 to the last integral in (jll|) . or rather to the upper and lower estimates 
that ([Tl~j) stands for, we see that 

y . / COsfA^'- 1 / 2 1 / l-CO S | 

G A (cos0,l)~~ l + -r ^ TTrrexp -c 



t J tA+i"^ * 

— (t + cos-l -^expl-c- 
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The last expression is what we need except for the first factor. However, since we consider 
A' > -1/2, 

t + cos-j >1, 0€[O,7r], 0<t<T, 



and, on the other hand, 



-A'-l/2 / n2 



t + cos-j <exp^e-J, 0€[O,tt], t > 0, 

for any e > 0, as can easily be seen by considering separately the cases 9 < it/ 2 and 9 > it/ 2. 
So this factor is insignificant, and Lemma 14.21 follows. □ 

Step 3. Let a, f3 > —1/2 be such that a + /3 is a dyadic number. Then A = Q + /3 + 1/2 > —1/2 
is also a dyadic number, so in view of Theorem 13.11 and Step 2 we may write 

G"^(cos 6, cos p) 



Cafi II G t/4 (usin-sin- + v cos - cos -, 1J dH a (u)dn.p(v) 
1 / „ , s / / 1 — u sin I sin ^ — u cos | cos 5 



<m Q (u) / ex P - c 2 — 2 - * — 2 - dn p ( v ) 



Applying now Lemma 14.31 as in Step 2, first to the integral in v and then to that in u, and 
observing that 

9 p 9 tp 2 9 — f -.2 

1 - sin - sin - - cos - cos - = 2 sin — ^— ~ (9 - tp) , 8, p € [0, vrj, 

we get 

(cos 0, cos </?) 

From this, ([9]) follows. In the next steps, we will remove the restriction that a + (3 is dyadic. 



Step 4. Suppose that A > is arbitrary. Then there exist e, e' > such that 2A — e and 2A + e' 
are dyadic numbers and A — e > 0. Applying Theorem 13.51 twice, with 8 = 0, j3 = A and either 
a = A — e or a = A, we obtain 

$ ej o(x)* e , (t/)G^(x,y) < ef( 2A+1 ~f)'Gf- e ' A (x,y), 

$ e ' ) o(^)$e' 1 o(y)G t A+e '' A (x,y) < e4( 2A + 1+ ^G A (*,y), 

for all x,y £ [—1, 1] and t > 0. This implies 

^ sin - sin ^0 G A+e ' A (cos 0, cos tp) < G A (cos 0, cos <£>) < ^ sin - sin ^ G A_e ' A (cos 0, cos ip), 

uniformly in 9, tp £ [0, tt] and < t < T; here and later on endpoint values are understood in a 
limiting sense, if necessary, and may be infinite. In an analogous way, we may vary the second 
type parameter and use Theorem 13.51 to get 

(cos-cos^ G t ,x+e (cos 9, cos tp) < G A (cos 9, cos p) < ( cos - cos ^\ G A,A e (cos 0, cos tp), 
uniformly in 0, £ [0, 7r] and < f < T. 
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Next we combine these estimates with those obtained in Step 3. For some positive constants 
ci and C2 and all 6,<p £ [0, ir] and < t < T, this leads to 

If 9 </j\-A-l/2/ (on-A-1/2 1 / (6 -iff 

7 (t + sin - sin — 1 1 + cos - cos — — = exp — ci 



(12) < G t A (cos6»,cos<y9) 

< F t (9,ip) e t + sin -sin- t + cos - cos - -= exp - c 2 



2 2/ V 2 2/ y/i V * 

where the auxiliary function F^ is defined by 

F t (0,<p) = min [ 1 + : * - , 1 + 



sm r; sin i| cos | cos 77 



It is easy to verify that 

F t (0 )¥ >)~l + — -. 



cos — 



Notice that the bounds in f)12|) coincide with those of Q, except for the factors involving Ft. 
In the next step, we will show how to deal with these factors. 

Step 5. We shall see how (fT2|) implies ([9]) with a = (3 = A > 0. Clearly, the factors F/ and 
1/F[' in (|12p are of importance only when 9 and (p are close to opposite endpoints of the interval 
[0,7r]. Indeed, we have 

27T 

F t (e,<p) e ' ' ~i~F t (e,<py, \e-<p\<—, o<t<r. 

Thus from now on we may assume that \9 — cp\ > 2n/3. Moreover, for symmetry reasons it is 
enough to consider the case 9 < it/ 3 and p > 2tt/3. 

Observe that under these assumptions, in the right-hand side of ([9]) only the exponential 
factor is significant, since it behaves like exp(— c/t) and the other factors are essentially contained 
between 1 and some negative power of t. So what we must prove is simply that 



(13) G A (cos 9, cos ip) ~~ exp 



for < t < T. We first verify f)13|) under the additional assumption that 9 > e c °/* or <p < 
7r — e _C0// *, where Co > is a sufficiently small constant. Clearly, either of these conditions 
implies cos ^ > e~ Co/t . Then for < t < T, 

F t (0, pf ~ ( 1 + —t— ) 6 < (1 + Te^Y < e c ^\ 
\ cos -^f- J 

and analogous bounds hold for Ft(6,(p) e '. This implies that if cq is taken small enough, the 
factors Fl and 1/F t e will be insignificant in (|12[) . and (|13p follows. 
Thus we proved the following. 

Lemma 4.4. Let A > 0. There exists a constant Co > such that ([9]) with a = /3 = A > holds 
for all 6, <p G [0, n] and < t < T, except possibly when 8 < e _C0// * and p > tt — e~ c °/* or vice 
versa. □ 

Finally, suppose that < e~ C0// * and > 7r — e~ C[,//t . By the semigroup property, 

/»7r 

G A t (cos 0, cos tp) = / G A (cos 0, cost/^G^cos '0, cos ip) (sin ■0) 2A+1 d0. 
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We split the interval of integration here into D\ = (0,e _c °/*), D2 = (e~ co ^ t ,ir — e~ c °/ 1 ) and 
D3 = (vr — e~ c °/*,7r), and we may assume that cq is small enough so that D2 D [ir /3,2ir/3]. 
Denote the resulting integrals by I\, I2 and Z3, respectively. To estimate I2, we apply Lemma 
14.41 getting 



r7, - e co/t r ' a ip\{ ip\ ( ip . (p\/ ip <p^- x ~ 1/2 



f'7T — e U' , 

' / ( 

Je-co/t [\ 



t + sin - sin — t + cos - cos — t + sin — sin — ][t + cos — cos — 
2 2/ V 2 2/ V 2 2/ V 2 2 



x - exp - c- — — — — (sin V>) + dip. 



t V * 

Here (6 — ip) 2 + {ip — ip) 2 ~ 1, which means that the exponential in the integrand makes all the 
other factors insignificant. Thus 

( c 

±2 ~~ exp I — - 

To bound I\, we apply the rough estimate of Theorem 13.71 to the two kernels in the integrand. 
We have 

c /i 



Zi<jT ° t- c (sin^) 2A+1 ^<exp(-|). 



Since the case of I3 is analogous, it follows that X\ +I2 and thus also G^(cos 8, cos 99) satisfy 
the estimate ([13]) . 

Altogether, we have proved the following. 

Lemma 4.5. The estimates © hold for a = (3 = A > and all 9, ip 6 [0, n] and < f < T . 

□ 



Step 6. From Lemma 14.51 with ip = 0, we get as in the last part of Step 2 
(14) G*(rafi 1 l)^4ap(-c?), #e[0,vr], < t < T, 



t A+i ^ V £ 

provided that A > 0; the case A = is covered by Step 1. Applying repeatedly Lemma 14.21 as in 
Step 2, we conclude that (114D holds for each A > —1/2, since the case A = —1/2 follows again 
from Step 1. 

Finally, we combine (|14p for arbitrary A = a + (3 + 1/2 > —1/2 with the reduction formula, 
as done in Step 3 for dyadic A. This establishes the short time bound of Theorem |A] for general 
a,f3 > -1/2. 

The proof of Theorem [A] is complete. 



5. The heat maximal operators 

We shall now consider d-dimensional Jacobi settings. Each of the semigroups T t a , 7^ 

and has a natural d-dimensional extension, see \18\ Section 2]. In particular, their kernels 

are simply tensor products of the corresponding one-dimensional heat kernels. Letting now 

a,/3 6 (— l,oo) rf denote type multi-parameters, we can use the same notation as before for the 

semigroups, their kernels and other related notions. The corresponding measure spaces will then 

be ([— 1, l] d , dg a ,p), ([0,Tv] d ,dfi a ^) and ([0, ir] d , d8), respectively, where 

d d 

dQa,/3 = (^) dQaifii, dji a ^ = (^) dflatfi., 
i=l i=l 

and d6 is the <i-dimensional Lebesgue measure in [0, 7r] d . 
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t>0 



and TT and T* with analogous definitions. Using Theorem we shall show that these 
operators satisfy weak type (1, 1) estimates in the corresponding measure spaces. 

Theorem 5.1. Let d > 1 and assume that a,j9 £ [— l/2,oo) d . Then 

(i) T* is bounded from L x {dg aj p) to weak L 1 (dg ct: i3); 

(ii) TT^ is bounded from L x (d[i a ^ to weak -L 1 (d// a a); 

(iii) T* is bounded from L 1 (d6) to weak L l (dO). 

An important consequence of Theorem 15.11 is the almost everywhere boundary convergence 
for the Jacobi semigroups applied to L l functions. Note that by the subordination principle, 
Theorem 15.11 implies weak type (1, 1) estimates, and thus also convergence results, for the multi- 
dimensional Poisson- Jacobi semigroups. 

We briefly discuss L p bounds for these operators. For a,/? £ (— 1,00)^, the boundedness of 
T,* in L p (dg a ^), 1 < p < 00, follows by Stein's general maximal theorem [221 Chapter 3]; see 
[TBI p. 346]. In the restricted range a,/3 G [— 1/2, oo)^, it can also be obtained by interpolation 
between Theorem 15, l( i) and the trivial boundedness in L°°. The case of 7** is much the same. 
In fact, 7** is controlled by T° , as can be seen from the proof given below, and so it inherits 
the LP mapping properties of T"'' 3 . The L p -boundedness of T* ,/3 is even simpler. Indeed, as 
pointed out in the proof below, Theorem lAl implies that when a,f3 G [— l/2,oo) d , the operator 
T*'^ is controlled by the standard maximal function in [0,7r] d and hence L p -bounded for p > 1. 
When a or (3 are not both in [— l/2,oo) d , the estimates of Theorem |A] can be expected to hold 
in the same form, and this suggests that the behavior of T* admits a similar anomaly to that 
occurring in certain Laguerre function settings and called pencil phenomenon |19j . 

Proof of Theorem \5.1\ From Theorem |A] it follows that for large t the three Jacobi heat kernels 
are bounded. Thus, from now on, we may consider only the maximal operators defined by taking 
suprema in the restricted range < t < 1. 

We first treat (iii) and (i) . In view of Theorem |A] and ([3|) , 



where c > depends only on a and /3. The right-hand side here is essentially the standard 
Gaussian kernel, and so T*'' 3 can be controlled by the Hardy-Littlewood maximal operator 
restricted to [0,7r] rf . Therefore, T°'^ is of weak type (1, 1). 

Next, we show that (i) follows from (ii). Observe that for / G L l (dg a ^) we have 



and consequently |T t /(cos#)| ~ \T" (/ o cos)(#)| for < t < 1. Thus (ii) implies (i). 

We pass to proving (ii). The kernel to be considered is then = ®^ =1 £7° , and we must 
estimate the supremum in < t < 1 of 






(15) 




9 G (0,vr) 
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with / G L l {dfi a ^). Here we may assume that / > 0. If we split each coordinate interval (0, 7r) 
in halves, the cube will be split in 2 d subcubes. It will be enough to consider in one of these, 
say (0, 7r/2) d , since the other subcubes can be treated like the one we select, if one places the 
Laguerre origin in the argument below at each corner of the large cube. 

Assuming thus G (0, Tr/2) d , we first do away with the integration over the set where 37r/4 < 
(pi < tt for some i. If 0, < it/2 and ipi > Sir /A, it follows from Theorem fAl that G^ 1 {9i, <Pi) < 1. 
For any nonempty subset A of {1, 2, . . . , d}, consider that part of the integral in (|15p taken over 
the set 

{p : 37r/4 < ifi < tt for i G A and < ipi < Stt/A for i £ A}. 

Carrying out the integrations over the interval (37r/4, tt) in each variable cpi with i G A, we 
can estimate this part of (|15p by 



n / S?" /,, (fl*,w)/A(^)dM^(vO. 



Here the product JJ/ is taken only over i G {1, . . . , d}\A, and (// consists only of the corresponding 
coordinates <pi. Further, f\{(p') is the result of integrating the given function / with respect to 
d^jL cti ^ i in (37r/4, tt) for each i G A. Finally, dfjf a g is the product of those measures dfi aii /3 i with 
i ^ A, each restricted to (0, 37r/4). In the extreme case A = {1, 2, ... , d}, the expression should 
be interpreted as Jr 37r u^d fdfia^, which certainly satisfies the weak type (1,1) estimate, so this 
case can be neglected. 

We thus arrive at an integral of the same type as that in (I15D but in lower dimension and 
with the important difference that the integration is now only over (0, 3-7r/4) in each variable. 
Summing up, this means that when estimating (1151) . we can restrict the integration to the cube 
(0,37r/4) d . What we must control is thus the operator 

Mf(9)= sup / 0€(O,7r/2) d , 

0<t<T7(0,37r/4) d 

for < / G L l (dn a fi). We will show that the weak type (1, 1) of M follows from the analogous 
property of a certain Laguerre maximal operator, which was proved by the authors in [IT] to 
satisfy the weak type (1, 1) estimate in the appropriate measure space. 
Assuming for a moment that d = 1, we observe that by Theorem El 

V ) ^ (t + ev)-^ 2 ^ exp ( - c^^) , 

uniformly in 6, ip G (0, 37r/4) and < t < 1. On the other hand, the one-dimensional heat kernel 
associated with the Laguerre system {xp% } considered in [TT] is expressed explicitly by 



K?(x, y) = j^- t exp ( - i coth t (x 2 + y 2 ) J (xy)- a I a 



xy 



2 sinh t 



for x,y,t > and a > — 1; here I a is the modified Bessel function of the first kind and order a. 
By means of the well-known asymptotics 



I a (z) ~ z a , z-^0 + , and I a (z) ~ z 1//2 exp(^), 
it is straightforward to check that 

K?(x, y) ^ (t + xy y a - 1 / 2 ^. exp ( - c t^Z. 



oo, 
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uniformly in x,y € (0, 37r/4) and < t < T, for any fixed T > 0. Thus we see that there exists 
a constant C > such that 

g?'P(e,<p)<K% t (e,<p), 9,<p G (0,3tt/4), 0<t<l. 

Moreover, the related Laguerre measure is drj a (x) = x 2a+1 dx and hence d[x a ^{9) ~ dr) a (9) for 
G (0,3tt/4). 

Coming back to arbitrary dimension d > 1 and taking into account the tensor product struc- 
tures of the Jacobi and Laguerre settings, we infer that for some C > 

g?'P(9,<p)<Kg t (6,<p), 9,<p£ (0,3^/4) d , 0<t<l. 

Further, the corresponding Jacobi and Laguerre measures are comparable on (0, 37r/4) rf . This 
reveals that M is controlled pointwise by the Laguerre maximal operator 

K?f(x) = sup / K?(x,y)\f(y)\ dr, a (y) 

t>0 J(0,oo) d 

restricted to the cube (0, 37r/4) rf . By [171 Theorem 1.1], satisfies the weak type (1, 1) estimate 
with respect to the measure space ((0, oo) d , dr) a ). We conclude that M is of weak type (1,1) 
with respect to ((0, 37r/4) rf , d/j, ai p), as desired. 

The proof of Theorem 15.11 is complete. □ 



Appendix: Poisson kernel estimates 

We complement the Jacobi heat kernel estimates by showing the following sharp bounds for 
the Poisson-Jacobi kernel in the Jacobi trigonometric polynomial setting. Clearly, this result can 
easily be transferred to the Jacobi trigonometric 'function' setting. In contrast to the preceding 
argument, the proof is based on an exact, positive representation of the kernel. 

Theorem 5.2. Assume that a,/3 > —1/2. Given any T > 0, we have 

H^(6, (p )^(t 2 + 6 2 + cp 2 ) ~ a ~ 1/2 (t 2 + (tt - 9f + (tt - tpf " '~ ' 2 



t2 + (0_p)2> 



uniformly in 9,ip G [0, tt] and < t < T, and 



K^,^)-exp(-t^±|±i 



uniformly in 9,<p G [0, tt] and t >T. 

The representation formula we shall use is |2Q, Proposition 4.1], which says that for a,/3 > 
-1/2 

(16) H t (9, V ) = c a , smh - J J (cQgh t __ i + ^ ^ M) v))a+fi+2 , 

with c Q)/ g = 2~ a ~ l3 ~ 1 / /j a ^(0, tt) and 

q(9, cp, u, v) = 1 — u sin - sin ^ — v cos - cos ^, 9, <p G [0, tt], u,ve[— 1,1]. 

This is based on the same product formula due to Dijksma and Koornwinder that we used in 
Section [3l The behavior of the double integral in (|16j) can be described by means of the following 
technical result. 
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Lemma 5.3. Let k > and 7 and v be such that 7 > v + 1/2 > 0. Then 

dn u (s) ^ 1 



/ hlil] (£> - £s) K (,4 - Bsp (D-B)*A"+i/ 2 (A- fl)7-«-i/2' < 5 < ,4 < 

Proof. We may assume that B = 1, since one can factor out a power of £? from both sides of 
the formula. The case f = —1/2 is trivial, so it is enough to consider i/ + l/2 > 0. Observe that 



.! (D - s) K (A - g)7 y (D-s) K (A-s)'r 7 (D - 1 + u) K (A - 1 + up 

Thus it suffices to analyze the last integral, which we denote by I. Now A > 1, and we consider 
the following two cases. 

Case 1: A > 2. Since in this case D — l + u~D — 1 and A — 1 + u ~ A — 1 ~ A for u G (0, 1), 
the conclusion is immediate. 
Case 2: 1 < A < 2. We split X as 




X1+X0. 



Then 



and 



(D - lp(A - 1)7 7 (£> - 1)«(A - 1)7-^-1/2 

1 „_1 /o ., . 1 



I„ < — — — _ / u "-l/2-7 <fo < 



(,D-1) K 7 A _ 1 ~ (D-1)*(A- l)7-"-V2 ' 

Since A ~ 1, this implies the conclusion. □ 

Proof of Theorem I5.£l We use (j!6p and apply Lemma 15.31 twice, first to the integral against 
cfflgtu), with the parameters v = j3, k = 0, 7 = a + /3 + 2, A = cosh | — usin|sin^, 5 = 
cos I cos ij, and then to the resulting integral against dH a (u), with the parameters v = a, 
k = /? + 1/2, 7 = a + 3/2, D = cosh ~, A = cosh | — cos | cos ^ and B = sin | sin This leads 
to the estimates 

sinh I 

(cosh I — cos I cos £) a+1 / 2 (cosh i _ s fn | s in £^ +1 / 2 (cosh I _ s i n I s i n £ _ cos I cos i£) 

1 

" ( cosh 1-1 + S in 2 ^ + sin 2 ?+£) Q+1 / 2 

1 sinh \ 



( cosh |_1 + sin 2 ^)-^-p) + sin 2 (*-»)+(*-?? ) W2 cosh \ - 1 + 2 sin 2 , 



As easily verified, the conclusion now follows. □ 
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